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RETARDATION OF PLATEAU-RAYLEIGH INSTABILITY: 
A DISTINGUISHING CHARACTERISTIC 
AMONG PERFECTLY WETTING FLUIDS 


JOHN MCCUAN 

Abstract. We consider a cylindrical film of fluid adhering to a rigid cylinder of 
fixed radius. The main result is to give the critical (maximum) length for which 
such a film of given thickness can be stable. The problem is considered both when 
the cylinder remains stationary and when the fluid and the cylinder are co-rotating 
at a fixed angular velocity. 


1. Young’s Equation 

The equation of Young asserts that the angle at which a fluid surface in equilib¬ 
rium meets homogeneous walls containing the fluid is constant (Figure |l](a)). This 



X 



Figure 1. Contact Angles 

assertion has been the subject of debate since it was introduced; Adamson [|, pg. 
360] writes concerning the microscopic contour of a meniscus, “The likely picture for 
the case of a nonwetting drop on a flat surface is shown in Figure 0(b). There is a 
region of negative curvature as the drop prohle joins the plane of the solid.” Pre¬ 
sumably, such a film (and perhaps the region of negative curvature) are not subject 
to variation. Nevertheless, if such a region exists even microscopically, then a math¬ 
ematical description of Figure 0](b) incorporating the macroscopic contact angle is of 
obvious interest. We offer such an explanation in an (admittedly) degenerate case 
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for which the film is subject to variation and (for this reason) does play a role on the 
macroscopic level. 


2. Perfectly Wetting Fluids 

Films of certain fluids assert their presence conspicuously due to their ability to 
transport the fluid. We have, here, helium at low temperature in mind, but there 
are other examples. From the point of view of equilibria, however, these fluids have 
typically been treated under the assumptions that the him is negligible and that 
Young’s condition prevails with 7 = 0 . Considering the typical thicknesses of such 
hlms (20-90 atomic diameters for helium) these assumptions may seem justihed. A 
recent experiment involving rotating nearly cylindrical hlms indicates otherwise. 
We now give a brief description of this experiment. 

In low gravity Q an annular container consisting of the region between two coaxial 
cylinders was partially hlled with a perfectly wetting huid. In this instance the inner 
cylinder was made of steel and the huid was silicone oil. Previous applications of the 
Young condition 7 = 0 and the assumptions described above predicted that under 
rotation about the axis of the cylinders (at certain speeds with certain volumes of 
huid) a stable axially symmetric equilibrium would be observed with cross-section 
resembling that shown in Figure ^ 



Figure 2. Axially Symmetric Doubly Connected Bubble 

What was in fact observed was no equilibrium at all, but the nearly cylindrical him 
on the inner cylinder wall thickened periodically and formed “pendant drops” which 
were hung outward, returning to the bulk of huid near the outer cylinder. In this 
way the huid was “pumped” by the capillary instability of the “him” on the inner 
cylinder in a cycle apparently against the centripetal force in part of its course. See 
Figure 3. 

It is this experiment which prompted the work to be presented below. 

^For this discussion we may assume zero gravity; in the actual experiment low gravity was 
simulated on the earth’s surface. 
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SHORT CYLINDRICAL CAPILLARY 


LONG CYLINDRICAL CAPILLARY 


SURFACE (STABLE) 


SURFACE (UNSTABLE) 


Figure 3. Thickening of Cylindrical Film 


3. A New Approach 


In the standard Young-Laplace theory of equilibrium capillary surfaces, for which 
the reader is referred to 0 , equilibrium surfaces S and their stability are determined 
by an energy functional of the form 



( 1 ) 


where a is the surface tension, |iS| is the area of the free surface, | W| is the area of the 
wetted region, X is the volume of fluid, and W is a potential depending on position; 
the Young equation is given by cosy = (3 where (3 is the “adhesion coefficient.” 
One might attempt to distinguish between various “zero contact angle” fluids in this 
framework by choosing an adhesion coefficient (3 > 1. It has recently been shown by 
B. White 0 that such an approach yields no difference (in equilibria nor stability) 
from the case (3 = 1. 

In the perfectly wetting case we suggest that the second term in (|^ be dropped 
altogether and a “Van der Waals” potential be incorporated in U. For containers 
with suitable geometry this additional potential depends only on the distance d to 
the container and, in any event, is required to become inhnite and negative as d \ 0. 
In this way, the issue of a contact line or a contact angle need never be addressed 
because all equilibria must wet every portion of the container. 

In the remainder of the paper we carry out an analysis of certain equilibria in a 
particularly simple container geometry. The results are rather surprising particularly 
in view of the historical introduction of the next section. 


4. Plateau-Rayleigh Instability 


Plateau in 1873 determined experimentally that the length at which a cylindrical 
column of liquid becomes unstable is between 3.13 and 3.18 times the diameter of the 


column.0 The theoretical value of vr for this multiple has been given by Rayleigh and 
more recently by Barbosa and do Carmo ii- We present a simple (and partial) 
derivation based on the second variation of (|l|). In this case | VV| =U = 0, and the hrst 


variation 6E gives rise to the equation of constant mean curvature 2H = A = constant 
to which the cylinder is a solution with H = l/2r. 


^This experiment was repeated in Q by G. Mason who obtained a value 3.143 ± .004. 





















4 


JOHN MCCUAN 


Calculating the second variation about the cylindrical solution one obtains 



where 9 and 2 ; are cylindrical coordinates on the cylinder S. Thus, using the criterion 
of Barbosa and DoCarmo p|, the cylinder will be unstable when L0 = — 

<pzz — has a non-positive eigenvalue whose eigenfunction satisfies 

/ 0 = 0. (3) 

J s 

This latter condition is the result of constrained volume. Solving L(j) = ficj) by a sepa¬ 
ration of variables 0 = A{6)B{z) with the boundary conditions A(0) = A{27r), A'(0) = 
A'{2ti)] 5(0) = 0 = B{1) one obtains eigenvalues 

— 1 

for fc = 0,1, 2,... and m = 1, 2, 3,.... 

I^k,m > 0 for A; > 0. Thus, instability can only arise for eigenvalues /xo,m- The 
smallest of these is /io,i = — 1/r^, but the corresponding eigenfunction 0o,i = 

sin(7r//)z does not satisfy and must be discarded. The next smallest is /io ,2 = 

— 1/r^ whose eigenfunction does satisfy and gives rise to an instability 
exactly when 

I > 27ir. ( 5 ) 

This is a symmetric variation, and at the length indicated in (^) the column of fluid 
will develop a neck which becomes smaller in radius until the column separates into 
two pieces. It is of note that in the experiment of G. Mason |^, the column of fluid 
was centered on a silica rod and the necking apparently proceeded until the column 
“ruptured” around the rod. According to our model such behavior is excluded for 
perfectly wetting fluids. 


5. Rotating Cylindrical Films 

We now carry over the analysis of the previous section to determine the critical 
length for the onset of instability of a cylindrical film of perfectly wetting fluid on a 
rotating, coaxial, rigid cylinder. 

The rotational potential energy is given by Ui = —{poi}‘^/2)r‘^, and we take a model 
Van der Waals potential approximating that for helium W 2 = —a/(B where d = t — Tq. 
A general form for Lcf) has been calculated by Wente and others : 

L(t) =-A,(t)-2{2H^ - K)4> + VU ■ N4> (6) 

where A is the intrinsic Laplacian on the free surface 5, H and K are the mean and 
Gauss curvatures of S and N is the outward normal to S. From this one obtains by 
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separation of variables as before 

dU — 1 

f^k,m = -^ + (") 

for = 0,1, 2,... and m = 1, 2, 3,.... The two smallest contributing eigenvalues are 

dU 1 47r2 


hog — 


dr P 

corresponding to the axisymmetric variation of the last section and 

_ dU TT^ 

+ F 


( 8 ) 


( 9 ) 


corresponding to a non-axisymmetric variation. The typical shape of dU/dr is shown 
in Figure 0 and evidently results in three cases. 


To 


a; = 0 



Figure 4. The Shape of dlA/dr 

VERY THIN FILMS: There is a value ri determined by the equation 

C) V l) 2 

or rf 


( 10 ) 


such that any cylindrical him (of arbitrary length) satisfying ro < r < ri is stable. 
For helium in the absence of rotation (|T0|) becomes 


3a 


(r — ro)"^ 


Thus, r = rcriticai satishes 


y/Scx 1 r r -— 7 

?^criticai = + -V4roV3a + 3a. 


(11) 


(12) 
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MEDIUM THICKNESS FILMS: Without rotation all cylindrical films with ri < r 
are unstable to the axially symmetric variation associated to /io ,2 whenever 

, , 27r 

I > lo,2 = 


dU 

dr 


In the presence of rotation there is a value r 2 > ri satisfying the equation 


(13) 


If ri < r < r 2 , then an axially symmetric instability sets in at / = /og. 

THICK FILMS (only in the presence of rotation): If r 2 < r, then the cylinder of 
radius r becomes unstable to the non-axisymmetric variation associated to /ii^i when 
I exceeds 

/ - "" 

“ 73' 

V dr 

NOTE: The only properties of the virtual (Van der Waals) potential U 2 which have 
been used are 

: (0) limr\,.o W 2 (r) = +cx). 

: (i) 

d^U2 


df2 


< 0 . 


: (ii) (When cu = 0) ^ decays faster than 1/r^ as r — 00 . 
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